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ABSTRACT 

In the standard scenario for spin evolution of isolated neutron stars, a young pulsar 
slows down with a surface magnetic field B that does not change. Thus the pulsar 
follows a constant B trajectory in the phase space of spin period and spin period time 
derivative. Such an evolution predicts a braking index n = 3 while the field is constant 
and n > 3 when the held decays. This contrasts with all nine observed values being 
n < 3. Here we consider a magnetic Held that is buried soon after birth and diffuses 
to the surface. We use a model of a growing surface magnetic Held to fit observations 
of the three pulsars with lowest n : PSR J0537—6910 with n = —1.5, PSR B0833—45 
(Vela) with n = 1.4, and PSR J1734—3333 with n = 0.9. By matching the age of 
each pulsar, we determine their magnetic field and spin period at birth and confirm 
the magnetar-strength field of PSR J1734—3333. Our results indicate that all three 
pulsars formed in a similar way to central compact objects (CCOs), with differences 
due to the amount of accreted mass. We suggest that magnetic field emergence may 
play a role in the distinctive glitch behaviour of low braking index pulsars, and we 
propose glitch behaviour and characteristic age as possible criteria in searches for CCO 
descendants. 

Key words: stars: magnetic field - stars: neutron - pulsars: general - pulsars: indi¬ 
vidual: PSR J0537—6910 pulsars: individual: PSR B0833—45 - pulsars: individual: 
PSR J1734—3333 


1 INTRODUCTION 

The magnetic field strength on the surface of neutron stars 
(NSs) spans a wide range: from B ~ 10 s —10 9 G for millisec¬ 
ond pulsars and NSs in low-mass X-ray binaries, through 
10 12 — 10 13 G for normal radio pulsars, to 10 14 — 10 15 G for 
magnetars. The primary method used to determine these 
magnetic fields is by measuring each pulsar’s spin period P 
and spin period time derivative P. Then assuming that the 
pulsar rotational energy decreases as a result of emission of 
magnetic dipole radiation, the surface field strength B at 
the magnetic polcQ is inferred, i.e., 

PP = 1b 2 47 B = 6.4 x 10 19 G (PP) 1/2 . (!) 

where 7 = An 2 R 6 sin 2 a/3c 3 1 = 4.884 x 

10 _4 ° s G -2 RqI jjd sin 2 a, R and I are the NS radius 
and moment of inertia, respectively, a is the angle between 


* Email: wynnho@slac.stanford.edu 

1 Note that a coefficient of 3.2 in eq. CD is often used in the 
literature, so that the inferred field in such a case is the field 
at the magnetic equator. Since we model field evolution at the 
magnetic pole, hereafter we only refer to the value at the pole. 


the stellar rotation and magn etic axes, and Re, = P/10 6 cm 
and J 45 = 7/ 10 45 g cm 2 (j Gunn fc Ostrikerl Il969l : see 

also ISpitkovskvl 120061 : IContopoulos et alJ 2014fl . Figure [T| 
shows the measured pulsar spin period and spin period 
derivative values taken f rom the ATNF Pulsar Catalogu<J3 
jManchester et ahl l200fit l. Also shown is a set of parallel 
lines which indicates the inferred magnetic field obtained 
from eq. ©• The other set of parallel lines indicates the 
pulsar spin-down or characteristic age t c = P/2P; r c is 
often used as a surrogate for the true age of a pulsar. 

The standard scenario for the rotational evolution of 
a pulsar is that it is born rapidly spinning (e.g., with an 
initial spin period Po in the millisecond regime) and rapidly 
slowing or spinning down, i.e., large P. This would place a 
newborn pulsar in the top-left region of Fig. |T] As it spins 
down, the pulsar moves for ~ 10 5 — 10 6 yr along a P-P 
path that tracks one of the short-dashed lines of constant 
magnetic field, evolving towards the bottom-right. This is 
because magnetic field diffusion and decay occurs on the 
Ohmic timescale 


2 http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Figure 1. Pulsar spin period P versus spin period time derivative P. Open squares denote pulsars whose values are taken from the 
ATNF Pulsar Catalogue; pulsars in a binary system and pulsars that have been observed to glitch are also noted. The nine stars denote 
pulsars with a measured braking index (see Table [IJ, with arrows indicating the direction each pulsar is evolving, as determined by their 
braking index. Closed and open circles and closed squares denote magnetars and isolated neutron stars (INSs) and central compact objects 
(CCOs), respectively. Short-dashed lines indicate characteristic age r c (= P/2P) and inferred magnetic field B [= 6.4 x 10 19 G (PP) 1 / 2 ]. 
Dotted line indicates the (theoretically uncertain) death line for pulsar radio emission. Solid lines indicate evolution trajectories (to age 
= 10 6 yr) where the magnetic field is buried initially and (from left to right) the initial field strength and spin period and burial density 
are [Po(G),Po(s),logpt>(g cm -3 )] = (3.4 x 10 12 ,0.015,11.0), (1.1 x 10 13 ,0.066,11.5), and (1.3 x 10 14 ,1.06,10.5), respectively. Long-dashed 
lines indicate trajectories where the magnetic field evolves according to eq. © and (from left to right) the initial field strength and field 
decay timescale are [Po(G),r(yr)] = (8 x 10 12 ,10 6 ), (3 x 10 13 ,10 6 ), and (3 x 10 14 ,10 5 ), respectively. 


4-7T cr c L 2 


TOhm = 


• 4 x 10 s yr ' 


V io 23 


1 km 


( 2 ) 


where a c is the electrical conductivity, L is the lengthscale 
over which decay occurs, and 1 km is the approximate size 
of the stellar crust (see Fig. EH); note that magnetic field 
changes can occur earlier for magnetars due to Hall effects 
which operate on a timescale rnaii = (4 nen e L 2 )/(cB ) ~ 
2 x 10 5 yr (p/10 12 g cnr 3 )(.B/10 14 G) _1 (L/1 km) 2 , where 
e is electron cha rge, n e is electron number den sity, and p 
is mass density jGoIdr eich fe Reiseneggen Il992l l (see also 
iGlampedakis et al.l 20111 1. In this work, we are primarily 
concerned with normal pulsars with B ~ 10 12 — 10 13 G and 
ages < a fewx 10 4 yr. Thus at times t <C rohm, the magnetic 
field does not change, and r c may be an adequate estimate 
of pulsar age (see Section [4]). However when t —> rohm, the 
magnetic field decreases, causing the efficiency of dipole ra¬ 
diation to decrease and P —¥ 0 [see eq. ©]• Simple examples 
of such evolutionary paths are shown by the long-dashed 


curves in Fig. [l] for different initial magnetic field Bo and 
decay timescale. In particular, we assume a magnetic field 
that evolves with time as 


m 


Bp 

1 + t/r’ 


(3) 


where r is the field decay timescale which can be taken to 
be approximately equal to rohm (or THaii for magnetars). 
Equation © mimics the results of numerical si mulations of 
magn etic field evolution in the crust (see, e.g., IColpi et al.l 
120001 . Observations and theoretical work seem to support 
the above scenario of an approximately constant magnetic 
field early in the li fe of a pulsar and a slowly decaying field at 
later times (e.g., IVigano et al.l 120131 find that the magnetic 
field is constant until an age of a few x 10® yr for Bp < 
10 14 G and ~ 10 s yr for Bo > 10 14 G, and llgoshev fe Popov! 
120141 find a field decay timescale of ~ 4 x 10 5 yr). 

On the other hand, there also exist observations that 
suggest that the magnetic field evolves, and especially fields 
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that grow, in young NSs. An important example comes from 
the (measured) braking index of pulsars. The second time 
derivative of the period P can be determined in a few pul¬ 
sars (where P is not dominated by timing noise), and P is 
conventionally expressed in terms of the braking index n, 
which is given by 


n 



(4) 


If pulsar spin-down is due to only magnetic dipole radiation 
and the field is constant, then eq. © yields a braking index 
n — 3. However, n < 3 for all pulsars with a measured P (see 
Table [T]). The low observed values of n can be attributed to 
a magnetic field that is increasing: allowing B to evolve in 
eq. 0 , one easily obtains 


n = 3 - 2~ = 3 - Tt-c^ = 3 - ~^P 2 , 
Bp B 7 B 3 


(5) 


where B is time derivative of B. For each pulsar with a mea¬ 
sured braking index, we denote its possibly evolving field 
in Fig. [T] by arrows direct ed along a pulsar’s trajectory in 
P—P phase space (see also lEspinoza et al.ll2011bl : lEspinozal 
l2013h . It is clear that about five of the nine pulsars are 
moving along trajectories almost parallel to (short-dashed) 
lines of constant B, and these are pulsars with braking in¬ 
dex between 2 and 3. Two pulsars (PSR J0537—6910 and 
J1734—3333) are clearly crossing lines of constant B, thus 
suggesting that their fiel ds are growing and J17 34—3333 is 
evolving into a magnetar (lEspinoza et alj |2011b"l. Note that 
we assume a constant a for simplicity [see eq. ©]• An evolv¬ 
ing a can produce similar spin evolution behaviour to one 
with B. However evidence for a varying a is uncertain (see, 
e.g..lLvne et al.i 2013, 2015h . and for example. iGuillon et al.l 
d2014ll find that either a is constant or the timescale for its 
variation is very long. 

In this work, we describe an alternative to the standard 
scenario for pulsar spin evolution described above, one that 
provides a physical mechanism for a grow ing magnetic field, 
using the model of 1 h 5 m r see also iMusli mov fe Pagel 


Il996r iGeppert et al. 19991 : Vigano &; Pond 20121 1. In brief, 
pulsars are born with a strong field, but this field was buried 
by an early episode of accretion and is slowly diffusing to 
the surface. In such cases, the surface field responsible for 
spin evolution [including n < 3 from eq. 0] is increasing 
at the current epoch. Since glitches and timing noise could 
be responsible for small changes in n (see Livm gstone et all 
l201lS : lAntonopoulou et ali 120151 : [Lvne et al.ll201.j i7 we only 
consider pulsars with a braking index less than half that 
predicted by the standard scenario (i.e., n < 1.5) as sources 
whose braking index re quires an explanation beyond the 
stand ard scenario (see lEspinoza et al.l 1201 lbl : iLvne et al.l 
120151 . and references therein, for discussion of othe r mod¬ 
els for low braking index; see also lHamil et al.ll2015l 'l. From 
Table[l]we see that this criterion is satisfied by three pulsars: 
PS R J0537—6910, B0833-45 (Vela), and J1734-3333; note 
that lMuslimov fc Pagel (|l996l l considered B0531+21 (Crab), 
B0540—69, and B1509—58 which have n > 2.1. In Section[2l 
we briefly describe the model for magnetic field evolution. 
In Section 02 we present our results and use the measured 
braking index and age of the three pulsars to determine the 
initial magnetic field and spin period of each pulsar and the 


amount of matter each accreted. In Section|4j we summarize 
our findings and discuss their implications. 


2 MAGNETIC FIELD EVOLUTION MODEL 


Our calculation of magne tic fi eld e volut ion follows that of 
lUrpin fe Muslimovl (1 19921 1 and iHol (l201lll . Here we provide 
a summary and describe updates (see morn . for details). 
The magnetic field is assumed to be buried deep beneath 
the surf ace by a post-su p ernova episode of hypercritical ac¬ 
cretio n dChevalierl 1 19891 : IGeppert et al.l 1 19991 : iBernal et al.l 
1201(1 l2013l l. The field then diffuses to the surface on a 
timescale that depends on burial depth. One might expect 
an anti-correlation between field growth rate and pulsar ve¬ 
locity since a smaller amount of mass will be accreted if the 
pulsar is moving at a greater velocity; observations t enta- 
tively support such a relation dGrinevdas fe Eksill2013l l. 

To determine the evolution of the buried magnetic field, 
we solve the induction equation 


dB 

dt 


= -V x 



V x B . 


( 6 ) 


Our interest is in the NS crust, which is predominantly in 
a solid state, and thus we neglect internal fluid motion. We 
take the surface field after NS formation, but prior to mass 
accretion, as the magnetic field strength at birth Bo. Accre¬ 
tion then buries and compresses this birth field . We assume a 
dipolar field in the stellar interior. In lHol (l201lh . we consider 
two field configurations: one in which the field is confined in 
the crust and another in which the field extends into the 
core. In the crust-confined case, the surface field grows at 
first but then decays at later times. I 11 the crust-core field 
case, the surface field strength grows until it saturates at 
the level of the core field, and the core field decays on a 
much longer timescale (> 10 6 yr, especially if the core is su¬ 
perconducting). Results for these two cases are qualitatively 
similar during the epoch of field growth. Here we consider 
only the latter case a nd hold the core field strength constant 
(see ILuq et akll2015l . for more results with a crust-confined 
field). A constant core field is justified since the Ohmic decay 
timescale [see eq. ©] in the core is longer than our times of 
interest (t •< 10 ® yr). 

We u pdated the magnetic field evolution code used in 
IHol (l201ll ! in primarily two ways . First, we use NS c rust and 
core models built with the APR (lAkmal et alll998lh BSk20, 
or BSk 21 (IPotekhin et al.ll2013l l models of the nuclear equa¬ 
tion of state (EOS). In the top panel of Fig. 0 we show the 
mass above a given density AM [= M — m(r), where M is 
total NS mass and m(r) is mass enclosed within radius r\. 
AM is an indication of the amount of accreted mass needed 
to bury the magnetic field to a given density. The bottom 
panel shows the crust depth (= R—r) as a function of density 
for different NS masses and EOSs. While AM(p) is weakly 
dependent on total NS mass and EOS, depth versus density 
is a strong function of M and EOS. For a given EOS model 
(e.g., APR), crust thickness decreases with increasing mass. 
Therefore to bury the magnetic field at a particular density, 
the field must be buried at a greater depth for a lower M. 
For different EOS models, we see that the result for BSk20 
is very similar to that for APR, and thus the field evolution 
timescale [which scales with depth, or L, as given by eq. (j2j 


© 2015 RAS, MNRAS 000, □□ 














































































4 Ho 


Table 1. Pulsars with a measured braking index n. Spin periods P and period derivatives P are taken from the ATNF Pulsar 


and from Ho et a 

■i <201511 

and references therein. References for age and braking index: 

ll iLvne et al. 19931. [21 Wane: & Gottheli 

199S 

: IChen et al. 

2006i, [3 

Middleditch et al. 20061. [41 Park et alj 20ld [51 Gradari et al. 

toll, [61 

Page et al.ll2009l: iTsuruta et al. 

200E 

. 171 Lvne et al. 199f 

. [81 Kumar et al.ll2012l. [9 

IWeltevrec 

e et al. 201 ll. [101 iGaensler et al. 

19991. [Ill Livine:stone &; Kaspi 

2011 

, [121 iHo & Andersson 20121 [131 Espinoza et al. 

2011b], [141 

Bocchino et al. 20051: ICamilo et al. 20061, [151 Rov et al. 20121, [161 

Blanton & Helfand 19961. [17] iLivineistone et al. 2007 

20li|. 


Pulsar 

SNR 

p 

P 

T c 

Age 

Braking 

No. of 

Typical 



00 

(s s- 1 ) 

(yr) 

(yr) 

index n 

glitches 

A9/Q 

B0531+21 

Crab 

0.0331 

4.23xl0 -13 

1240 

961 

2.51(1) [1] 

25 

10“ 9 - 10~ 8 

J0537—6910 

N157B 

0.0161 

5.18xl0 -14 

4930 

2000[t 3 qoo [2] 

-1.5(1) [3] 

45 

10- 7 

B0540-69 

0540-69.3 

0.0505 

4.79xl0 -13 

1670 

1000t®®° [4] 

2.087(7) [5] 

1 

10~ 9 

B0833-45 

Vela 

0.0893 

1.25xl0 -13 

11300 

11000^°° [6] 

1.4(2) [7] 

19 

10-® 

,11119-6127 

G292.2—0.5 

0.408 

4.02xl0 -12 

1610 

7100t 3 °° 0 [8] 

2.684(2) [9] 

3 

10" 9 - 10" 6 

B1509-58 

G320.4—1.2 

0.151 

1.53xl0 -12 

1570 

< 21000 [10] 

2.832(3) [11] 

0 

— 

.11734-3333 

G354.8—0.8 

1.17 

2.28xl0 -12 

8130 

> 1300 [12] 

0.9(2) [13] 

1 

10- 7 

J1833-1034 

G21.5—0.9 

0.0619 

2.02xl0 -13 

4850 

10001”° M 

1.8569(6) [15] 

4 

10- 9 

.11846-0258 

Rest even 75 

0.327 

7.11xl0 -12 

728 

1000l 3 gg° [16] 

2.65(1) [17] 

2 

10" 9 - 10-® 


for the Ohmic timescale] for BSk20 and APR will be compa¬ 
rable. Depth at a given density for NSs built using BSk21 is 
larger than that for NSs built using APR, and thus the field 
evolution timescale for BSk21 will be longer than that for 
APR. Note that, in addition to depth variations between 
the three EOS models considered here, crust composition 
for each model is different, which in turn produces differ¬ 
ent elect rical conductivitie s . We ment ion that the previou s 
works of iMuslimov fe Page! dl996li and lGeppert et al.l (Il999l ) 
consider different EOS models than those studied here. 

The second update is that we use CONDUCT110, which 
implements the latest advancements in calculating electrical 
conductivities (IPotekhin et al.l 120151') . We assume no con¬ 
tribution due to impurity scattering since this only be¬ 
comes important at high densities and low temperatures. 
We checked that there are no noticeable changes for a uni¬ 
form impurity parameter Q i mp < 1 (where a c oc Q[^ p ), 
which is the relevant regime for the crust of isolated NSs, 
in contrast to that of NSs accreting from a binary compan¬ 
ion. Recent works examine t he effects of larger Qj mp on spin 
and magnetic field evolution llPons et al.ll201 j : IVigano et al.l 
1 201 .'ll : iHorowitz et al.l |2015h . However, this occurs due to 
pasta phases near the crust-core boundary at densities ~ 
10 14 g cm ~ 3 , and its effects only become im portant after 
~ 10 s yr dPons et al.ll2fll3i : IVigano et al.ll2013l ). 

3 FIT TO LOW BRAKING INDEX PULSARS 

The age of the three pulsars (PSR J0537—6910, B0833—45, 
and J1734—3333) with n < 1.5 is in the range of sa 1 T6 kyr 
(see Table [lj. Therefore, we seek a magnetic field growth 
timescale on this order. From eq. © and Fig. [2] we esti¬ 
mate that the magnetic held should be buried at a den¬ 
sity ~ 10 11 g cm -3 , which corresponds to an accreted mass 
AM ~ HT S M 0 . 

Beginning with an initial condition which has the mag¬ 
netic field Bq buried at density pb, the induction equation 
[eq. ©] is solved to obtain the radial profile of magnetic 

3 http://www.ioffe.ru/astro/conduct/ 



log p (g cm 3 ) 


Figure 2. Accreted mass AM (top panel) and depth (bot¬ 
tom panel) as a function of density p. Vertical dotted lines in¬ 
dicate boundaries between envelope and crust (taken to be at 
10 10 g cm -3 ), outer and inner crust (at Pdrip ~ 4x 10 11 g cm -3 ), 
and crust and core (at baryon density n = 0.09 fm —3 ). 


field as a function of time. Integrating eq. ©, we obtain the 
evolution of the spin period 


P(t) = 


Po+7/ B(R,t') 2 dt' 
Jt 0 


-I 1/2 


(7) 


where Po is the initial pulsar spin period and B(R,t) is the 
magnetic held at the NS surface at time t and is given by our 
solution to the induction equation. From P(t), we calculate 
its hrst and second time derivatives P(t) and P(t), respec¬ 
tively, and corresponding braking index n given by eq. ©. 
For each of the three pulsars, we vary the three initial con- 
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ditions [Bo(G),Po(s),logpb(g cm~ 3 )] until the resulting P, 
P, and n match those of the pulsar at its current age. It is 
worth pointing out (see also lHo|[2013b f that fitting to ob¬ 
servables P and P is equivalent to fitting to r c , B as given by 
eq. Q. or spin-down lu minosity since these are all formed 
from P and P (see, e.g.. iMuslimov fc Pagelfl996l 'l. while n is 
an independent parameter since it includes P. 

Figure [T] shows results for a set of initial conditions 
that fit each of the three pulsars: (Bo,Po,log pb) = (3.4 x 
10 12 ,0.015,11.0) for J0537—6910, (1.1 x 10 13 ,0.066,11.5) for 
B0833-45, and (1.3 x 10 14 ,1.06,10.5) for J1734-3333. Each 
trajectory is labelled with the approximate age at which the 
calculated P, P, and n match their corresponding observed 
values. For J0537—6910, we obtain an age of 1.95 kyr, com¬ 
pared to its known age of 2^ 3 kyr. Because of its young age, 
the current spin period P = 16 ms is not much different from 
the initial spin period Po = 15 ms. For B0833—45, we obtain 
an age of 10.2 kyr, compared to its known age of lljlij 6 kyr. 
The spin period in this pulsar has increased significantly 
from 66 ms to the current 89 ms due to its much older age. 
For J1734—3333, we obtain an age of 2.07 kyr, compared 
to its minimum age of > 1.3 kyr. Despite its similar age to 
J0537—6910, J1734—3333 underwent noticeable spin-down 
due to its much stronger magnetic field. Since B ~ 10 14 G 
for J1734—3333, Hall effects and anisotropic conductivities 
may need to be taken into account, although these should 
not change our conclusions at a qualitative level . We also 
note the work of lGourgouliatos fc Cummingl d2015l l, who use 
numerical simulations to model Hall drift effects and obtain 
field growth with multipolar toroidal fields of the order of 
10 14 G that can produce braking indices matching those ob¬ 
served. 

Figure [3] shows relationships between the three initial 
conditions, Bo, Po, and pb, as well as how initial values com¬ 
pare to present values. For a pulsar at a given age, the in¬ 
ferred density at which the magnetic field is initially buried 
depends on NS mass and nuclear EOS model, while Bo and 
Po do not depend significantly on M or EOS. The thicker the 
crust (see Fig. 0 , the shallower the field needs to be buried 
to match current values of P, P, and n. The shallower the 
field is buried, the closer the current magnetic field and spin 
period are to their initial values. Note that B is evolving 
towards Bo, while P is evolving away from Po. 

More accurate age determinations would lead to tighter 
limits on the burial density (and hence mass accreted, via 
Fig. 0 and initial magnetic field and spin period. Cur¬ 
rent age estimates allow us to constrain the magnetic held 
strength at birth to within a factor of about three, i.e., Bo « 
(3 - 7) x 10 12 G for J0537-6910, B 0 « (0.9 - 1.3) x 10 13 G 
for B0833-45, and B 0 ~ (1 - 3) x 10 14 G for J1734-3333. 
These field strengths fall within the lognormal distributions 
determined from population synthesis studies, which find 
an average and width cr of log Bp = 12.95 ± 0.55 in the 
case where there is no field decay dFaucher-Giguere fc Kaspil 
120061 : IGuillon et alJ l2014li and log Bo = 13.25 ± 0.6 in 
the case of where there is (model-de pendent) field decay 
lIPopov et al.ll20i(il : IGuillon et al.ll2014T ). The initial spin pe¬ 
riods are somewhat shorter for J0537—6910 and B0833—45 
and much longer for J1734—3333 than the average found in 
these population synthesis studi es, i.e., Po = 0.30 ± 0.15 s 
( Faucher-Gi guere fc Kaspil l2006ll and Pp = 0.25 ± 0.1 s 
( Popov et al. 20101 1. However Guillon et alJ (1 201 ill find a 


much broader distribution, such that even Po as 1 s for 
J1734—3333 is within 2<r of the average value (see also 
llgoshev fc Pot)O\1l20i ill ). 


4 DISCUSSION 

The standard theoretical scenario for spin evolution of iso¬ 
lated pulsars is one in which a pulsar loses its rotational 
energy via dipole radiation and slows down over time with 
a constant surface magnetic field, and this scenario works 
well to explain pulsars at an early age (< 10 s yr) and 
pulsars with braking index n as 3 (see Table |T]). Here we 
study an alternative scenario, one in which the intrinsic 
magnetic field of a pulsar is buried by accretion soon af¬ 
ter its birth in a supernova. Our scenario works alongside 
the standard scenario, by providing a natural explanation 
for pulsars with n < 3. We fit observations of the three 
pulsars whose braking index n < 1.5, i.e., PSR J0537—6910, 
B0833—45, and J1734—3333. We find that the mass required 
to bury the magnetic field is AM ~ 10 _5 Mq. Another 
requirement is that the timescale over which this mass is 
accreted must be sh orter than the field diffusion timescale 
iGeppert et al.l[l999l i . given approximately by eq. (0. Thus 
most newborn NSs probably fall within the standard sce¬ 
nario, with very little or no accretion. A relatively few, like 
the three examined here, quickly accreted enough mass to 
bury their magnetic field. Generally, the range of accreted 
mass could be quite large. Thus a third formation channel 
is one in which a large amount of matter is accreted. Such 
is possibly the case f or the central compact objects (CCOs) 
studied in iHol (1201 lh . where AM > 10 ~ 4 Mq. CCOs have 
age < a few x 1 0 4 yr, and their surface magnetic field is 
~ 10 10 — 10 11 G (lHalpern fc Gotthelfl l2oiol : iGotthelf et all 
l2013al:lHoll2013al! . 

Unification of different observational classes of NSs, 
such as magnetars, CCOs, and normal radio pulsars, via 
evolution of a NS from one class to another (iKas aiiia 
IPqpov et alJl2C)i(ih . would help allevia t e the NS-supernova 
birthrate problem (iKeane fc Kramer! l2008l l. One pulsar 
studied here that is of particular interest in this regard is 
J1734—3333. Based on its current inferred magnetic field 
and braking index, this pulsar appears to be moving in 
Fig. 0 from the region p opulated b y normal ra dio pulsars 
into that of magnetars (lEspinoza et all 1201 lbl l. We find 
that it has a magnetar-strength magnetic field at birth, i.e., 
Bo ~ (1 — 3) x 10 14 G, with the exact value dependent on the 
age of the pulsar (see right-hand panel of Fig. 0 . In the case 
of shallow field burial and no significant field decay during its 
life thus far, the magnetic field of J1734—3333 is nearly at its 
birth value and will reach it at an age as 10 4 yr. We point out 
that our projection of the trajectory of J1734—3333 is one 
that involves a brakin g index which changes over time (cf. 
lEspinoza et all 1201 lbli . In particular, a detectable change 
from n = 0.9 to 1, based on the current uncertainty of 0.1, 
takes about 70 yr. J0537—6910 is perhaps a more promising 
target for observing a braking index change: a change from 
n = —1.5 to —1.4 takes about 20 yr. For completeness, the 
braking index change for Vela from n = 1.4 to 1.5 takes 
about 400-500 yr. The difficulty lies with the much larger 
braking index uncertainty for these three (low braking in¬ 
dex) pulsars compared to that of other pulsars (see Table |T|). 
For example, evolution of braking index and third time 
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Figure 3. Initial magnetic field Bo (top), initial spin period Pq (centre), and burial density p\, (bottom) as a function of age for 
PSR J0537—6910 (left-hand panel), B0833—45 (centre panel), and J1734—3333 (right-hand panel). Horizontal dotted lines denote the 
current magnetic field B and spin period P of each pulsar, where the former is given by B = 6.4 X 10 19 G (PP) 1 / 2 . Vertical solid and 
dotted lines indicate the nominal age and age range, respectively, of each pulsar (see Table [TJ , except in the case of J1734—3333, the 
vertical dotted lines indicate the age lower limit. 


derivative of spin period is observed for Crab ( Lvne et al.l 


120151 ) and B1509—58 (jLivingstone fe Kaspillioil ). 


It is possibly noteworthy that two (J0537—6910 and 
Vela) of the three pulsars with n < 1.5 undergo regular, 
large-amplitude glitches in their timing behaviour. The third 
(J1734—3333) recently had a large glitch, as reported in 
the Glitch Catalogu<0 (lEspinoza et al] l2011al ). The other 
pulsars with 1.8 < n < 3 either are not seen to glitch 
(B1509—58), have small amplitude glitches (e.g., Crab), or 
have glitches whose amplitude varies greatly (e.g., the mag- 
netar J1846—0258) (see Table [T]). The distinctiveness of 
J0537— 6910 and Vela glitches is discussed bv lEspinoza et all 
ll2011al) . who also show that glitc h size has a bimodal dis¬ 
tribution (see also lYu et al.ll2013l ). The possible connection 
between low braking inde x and regu l ar, la rge glitches in cer¬ 
tain pulsars is noted by lEspinozal (120131 ) , who report an¬ 
other three pulsars (with n = —1,2, and 2) that could 
be in this group. These large amplitude spin-up glitches 
are of great interest since they may be revealing proper¬ 
ties of the neutron sup e rfluid in the NS llBavm et al. ] [ 19691: 
lAnderson fe Itohl Il975l : lAlpar et al.l 1 19841 ). The regularity 
of these similarly-sized glitches is thought to be the re¬ 
sult of the pulsar tapping and exhausting the entire an- 


ner crust (Link et al. 1999: lA rdersson et al. 

20121: Charnel 

2012 

: Piekarewicz et al.l 20141: Ho et al.l 2015 

: Hooker et al.l 

2015 

: Steiner et al. 20151). Our simulations of magnetic field 


diffusion from the (inner and outer) crust to the surface sug¬ 
gest that perhaps this motion could be involved in triggering 
glitches of the type seen in low braking index pulsars. 

If regular, large glitches are a symptom of a previ¬ 
ously buried magnetic field, then glitch activity could be 
used as a criterion for searches for descendants of CCOs. 
Previous searches for CCOs and their descendants focus 
on the region in P-P phase space where known CCOs re¬ 
side (see Fig. [lj. Th ese searches have thus far not found 
definitive candidates dCotthelf et al.ll2013bl: Bogdanov et al 


120141 : ILuo et al.l [2015!) . As we show here (see also lLuo et al 


4 http://www.jb.man.ac.uk/pulsar/glitches.html 


120151) . pulsars with an emerging magnetic field move rapidly 
through this region of P-P, and thus the likelihood of dis¬ 
covery here is l ow. This can also explain the relative p aucity 
of pulsars here (jHalpern fc Cotthel l20lrt : lKaspil201(il ). Once 
their intrinsic fields reach the surface and these pulsars 
evolve to join the majority of the NS population, they might 
still be distinguished by their glitch activi ty. After all, glitch 
size and activity peak at a ge ~ 10 4 yr llMcKenna fe Lvna 
1 1990l : [Espinoza et al.ll201 1a ), and glitches only occur in pul¬ 
sars with r c < 10 Myr ( Espinoza et al.ll2011al ). 

Finally, it is well known that pulsar characteristic age 
r c is often discrepant with true age, and thus the former 
can be an unreliable estimate of the latter (e.g., the case 
with CCOs). For seven of the nine braking index pulsars, 
the characteristic age is within a factor of about two of the 
true age (see Table Q]), although some of the true age deter¬ 
minations are likely biased towards r c . Thus t c is a relatively 
good age estimate for these types of sources and when these 
pulsars are near maximum P and post-maximum. Thus low 
r c could be another possible criterion in searches for CCO 
descendants. 
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